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Abstract 

We attempted to control the incorporation of twin boundaries in self-catalyzed GaAs nanowires (NWs). 
Self-catalyzed GaAs NWs were grown on a Si substrate under various arsenic pressures using molecular beam 
epitaxy and the vapor-liquid-solid method. When the arsenic flux is low, wurtzite structures are dominant in the 
GaAs NWs. On the other hand, zinc blende structures become dominant as the arsenic flux rises. We discussed this 
phenomenon on the basis of thermodynamics and examined the probability of twin-boundary formation in detail. 
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Background 

III-V compound semiconductor nanowires (NWs) have 
been attracting significant attention as fundamental 
structures of novel optical and electronic devices. Es- 
pecially, vapor-liquid-solid (VLS) NWs grown on Si 
substrates have been investigated for optoelectronic 
integrated circuits [1,2]. 

Previously, we succeeded in growing self-catalyzed 
GaAs NWs on Si substrates using molecular beam epi- 
taxy (MBE)-VLS method which is a combination of 
MBE and VLS method [3]. However, twin boundaries 
formed in the NWs during growth. The occurrence of 
twin boundaries is controlled by various methods such 
as control of supersaturation, growth temperature, and 
diameter [4-8]. To control the twin boundaries, we fo- 
cused on the pressure of arsenic because self-catalyzed 
GaAs NW growth was nearly independent of Ga pres- 
sure, and arsenic flux plays an important role in the 
growth mechanism [3]. In addition, arsenic solubility in 
Ga solution is very low [9]. This means that the degree 
of supersaturation depends on arsenic pressure only. 



Methods 

Self- catalyzed GaAs NWs were grown on a (lll)Si sub- 
strate by MBE-VLS method. The growth temperature 
was 580°C. When the arsenic flux varied from 5.0x10 r 6 
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to 1.9x10" Torr, the diameter and the length of the 
obtained GaAs NWs varied from 90 to 30 nm and 0.5 to 
3.5 (im, respectively. We have previously discussed the 
tendency of the diameter and length to depend on the 
arsenic flux [3]. The obtained NWs were observed by 
transmission electron microscopy (TEM). 

Results and discussion 

TEM images of GaAs NWs are shown in Figure 1. In 
the case of low arsenic flux (5.0 x 10" 6 Torr), wurtzite 
(WZ) structures are dominant as shown in Figure la. 
On the other hand, when the arsenic flux is high, 
Figure lc shows that segments between the twin bound- 
aries become large, and zinc blende (ZB) structures are 
dominant. The segment size follows the time between 
successive twin-crystal nucleation events. Furthermore, 
distribution of the segment sizes is exponential from a 
stochastic point of view [10]. 

Figure 2 shows the distribution histogram of the seg- 
ment sizes obtained from Figure 1. We assume a segment 
size of x and fitted the histogram with an exponential 
curve of the form exp (~x/a) to estimate the expectation 
value of segment a. The estimated expectation values were 
1.1 and 6.5 monolayers at arsenic fluxes of 7.0 x 10" 6 and 
1.9 x 10" 5 Torr, respectively. The reciprocal of the expect- 
ation value a" 1 is equivalent to the probability of occur- 
rence of twin-crystal nucleation. The probabilities were 
approximately 90% (7.0 x 10" 6 Torr) and 15% (1.9 x 10" 5 
Torr). The dependence of the probability on the arsenic 
flux is in good agreement with [11] and [12]. 
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Figure 1 TEM images of NWs grown under various arsenic fluxes, (a) 5.0 x 10 6 , (b) 7.0 x 10 6 , and (c) 1.9 x 10 5 Torr. 



To understand these phenomena, we calculate the de- 
gree of supersaturation and estimate the probability of 
twin-boundary formation, following the procedure pre- 
sented by Glas [13]. The degree of supersaturation A\i is 
as follows: 



A{JL 



■2# 



GaAs 



(i) 



where ^ Ga As is the half chemical potential of GaAs crys- 
tal nucleation; ^a an d ^a s are the chemical potentials of 
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Figure 2 Histogram and distribution of segment size obtained 
from Figure 1. Arsenic fluxes are (a) 7.0 x 10" 6 and (b) 1.9 x 10 

~ 5 Torr. Solid lines are exponential fitting curves. 



liquid gallium and liquid arsenic, respectively. We as- 
sume that the arsenic adatom was liquid arsenic existing 
on the gallium droplet. However, ^a an d f*As include an 
interaction between gallium and arsenic in pure gallium 
and arsenic chemical potentials an d f*As> respect- 
ively. Therefore, ft x (X = Ga or As) is given by 



+RTln(a x ) 



(2) 



where R is the gas constant, T is the growth 
temperature, and is the activity of X in the liquid 
phase. At T 0 = 298.15 K, when we adopt an enthalpy per 
mol in the X, the solid phase of h x % the degree of super- 
saturation is given by 

Afi = RTln(a L Ga ) +RTln(a L As ) + (fa - h G %) 

+ (^s-Co)-2(^GaAs " 0.5< 0 - 0.5< 0 ) . 

All terms except the enthalpies depend on the growth 
temperature. The third, fourth, and fifth terms are the 
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differences between the chemical potential and enthalpy 
of the pure gallium droplet, arsenic droplet, and GaAs 
crystal Their values are tabulated in the paper of Ansara 
et al [14]. The first and second terms are concerned 
with interactions in the alloy. They are denoted by the 
atomic concentration c x and the interaction parameter 
for gallium and arsenic o) G ^ As as follows: 

*rin(4a(As)) = ^ r ^( C Ga(As))+4s(Ga) W Ga,As 



:2?r In (c Ga (As)) + 

X (§Ga - §As) 



As(Ga) V GaVAs 
CGaVba + CAsVXs 
1.256 Xl0 5 ( XGa 



XAs 



where VG a(As ), Xg*(As)> and £ Ga(As) are the molar volume 
in the liquid phase, Pauling electronegativity, and Hil- 
debrand solubility parameter, respectively. From these 
equations, we can obtain the relational expression for 
the supersaturation dependence on growth temperature 
and arsenic concentration in the gallium droplet. 

Figure 3 shows the supersaturation as a function of the 
arsenic concentration in the gallium droplet at 580°C. In 
this calculation, the supersaturation per atom Ap/N A , 
where N A is Avogadro's number, is adopted. In Figure 3, 
the equilibrium arsenic concentration is approximately 
0.06% at Ap/N A = 0. This value is in fair agreement with 
the arsenic concentration of approximately 0.1% in the 
gallium droplet at 580°C [9]. Therefore, when the arsenic 
concentration is positive because of an oversupply of 
arsenic atoms greater than 0.06%, GaAs VLS growth 
may occur. 

From supersaturation A(aIN Ai we can obtain the 
probability of twin-boundary formation. When arsenic 
adatoms on the gallium droplet surface diffuse to the 
three-phase boundary of vapor, liquid, and solid phases, 
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Figure 3 Supersaturation as a function of arsenic concentration 
in gallium droplet at 580°C. 



we assumed that the GaAs crystal nucleus would be 
formed into a rhombus shape with side length r at the 
vertex of the nanowire top surface. Figure 4 shows the 
growth model. When the ZB and WZ structures form, 
the amount of Gibbs free energy change in this growth 
system is given by 



AG 



ZB(WZ) 



2D, N A 

+ 2rh |xsl,zb(wz) + £Vsv,zb(wz) _ Ylv sm /? j 

+ — s-Xsn 



(4) 



where H is the solid volume of the Ga-As pair (4.51 x 
(3) 10" 29 m 3 ), h is the height of crystal nucleus, ft is the 
contact angle of the droplet, y LV (y SN ) is the interface en- 
ergy between the droplet and the vapor phase (between 
the crystal nucleus and the NW top surface), and ^sl,zb 
( WZ) and Xsv,zb(wz) are the interface energies of the GaAs 
ZB (WZ) crystal nucleus top and side surfaces. Obviously, 
h is the (111) GaAs lattice spacing, which is 0.32639 nm. 
Xsv,zb and Xsv,wz are the surface energies of (l-lO)GaAs 
and (ll-20)GaAs, which are 0.62 and 0.54 J/m 2 , respect- 
ively [15]. The interface energy y LV depends on 
temperature, and we use the following relation: y LV = 0.708 
- 0.66 x 10" 4 x (T- 303) (J/m 2 ) [16]. If the crystal nucleus 
is ZB, y SN = 0. On the other hand, in the case of the WZ 
crystal nucleus, y s ^ is 0.023 (J/m 2 ), which is half of the 
GaAs stacking fault energy [17]. We assumed that the 
contact angle ft is 45° [18]. Since the NW side surface has 
a certain asperity, we adopted the parameter C (0 < C < !)♦ 
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Figure 4 Growth model adopted in the calculation. 
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Therefore, the maximum AG Z b(wz)* and the critical nu- 
cleus of Equation 4 for r are obtained as follows: 



AG 



ZB(WZ) 



'ZB(WZ) 



|ksl,zb(wz) + C Ksv,zb(wz) Yin sm /?| 

2 ^{>"SL,ZB(WZ) + CXsv,zb(wz) ~~ KLV sin ^} 



(5) 



At steady state, it is known that the probability of crys- 
tal nucleation /zb(wz) is proportional to the S ZB (wz) re- 
gion, where a crystal nucleus can form, and Zeldovich 
factor Z ZB(WZ ) [19]. The relation is given by 

(^^ZB(WZ) \ 
kf / ^ 

where h is the Boltzmann constant, Szb(wz) = 
3 y/3 ^zb(wz)\ > anc * tne Zeldovich factor is 



z Z B(wz) = {v^(^-rs N )/^r} 



1/2 



Therefore, the probability of twin-crystal nucleation P is 

r, IwZ 



JzB + /WZ 



ZwzSwz exp(-^p) 



Z zb S Z b exp (- + Z wz «S W z exp (- ^|p) 



(7) 



This equation indicates that the probability depends 
on the growth temperature and arsenic concentration. 
When we assume that ( is 0.8 at 580°C, the relation be- 
tween the probability P and arsenic concentration c As is 
shown in Figure 5. As the droplet arsenic concentration 
increases with increasing arsenic flux, the probability of 
twin-crystal nucleation decreases. This phenomenon 
agrees with the experimental results. 

By using the probability of twin-crystal nucleation in 
Figure 2, we can calculate the arsenic concentration and 
supersaturation per atom. When the arsenic fluxes are 
7.0 x 1(T 6 and 1.9 x 10 5 Torr, the arsenic concentra- 
tions are 0.11% and 0.38%, and the supersaturations per 
atom are 53 and 143 meV, respectively. These supersatur- 
ation values are smaller than those of Au-catalyzed GaAs 
NWs (230 to 1,570 meV) [20]. This difference might be 
due to the difference in the side facet surface. Glas et al. 
adopted the {111} and {1-100} facets in the ZB and WZ 
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Figure 5 Relation between the probability of the twin-crystal 
nucleus and arsenic concentration at 580°C (C = 0.8). 



structures in their calculations, respectively. In addition, 
there might be a difference in the diffusion length be- 
tween gold and gallium droplets. From the obtained ar- 
senic concentration, we estimate the critical nucleus. 
When the arsenic fluxes are 7.0 x 1(T 6 and 1.9 x 10" 5 
Torr, the critical nuclei of ZB (WZ) are 1.1 (0.3) and 0.4 
(0.1) nm, respectively. This means that the increase of 
arsenic flux decreases the critical nucleus and increases 
the growth rate. In the case of high arsenic flux, the size 
difference between critical ZB and WZ nuclei is small 
compared with the case of low arsenic flux. This means 
that the ZB structure appears easily as the arsenic flux 
increases. Therefore, we could improve the comprehen- 
sion of the growth mechanism in the self- catalyzed GaAs 
NWs. This comprehension might support a techno- 
logical feasibility of a novel device like twin-plane ID 
superlattices [21]. 

Conclusions 

Self-catalyzed GaAs NWs were grown on a (111) Si sub- 
strate by MBE-VLS method under various arsenic fluxes. 
From the TEM observations, we found that the segment 
size between the twin boundaries depends on the arsenic 
flux. In order to understand this phenomenon, we 
attempted to calculate the degree of supersaturation and 
estimate the probability of twin-boundary formation. 
When the supersaturation increased with increasing ar- 
senic flux, the size difference between the critical ZB 
and WZ nuclei decreased. As a result, the ZB structures 
were easier to obtain as the arsenic flux increased. This 
qualitatively explained the experimental results and the 
high probability of the incorporation of twin boundaries. 
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